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Abstract. 

^ Several quantities related to the Zernike circle polynomials admit an ex- 

t^^ pression as an infinite integral involving the product of two or three Bessel 

^ functions. In this paper these integrals are identified and evaluated explic- 

(^ itly for the cases of (a) the expansion coefficients of scaled-and-shifted circle 

1^ polynomials, (b) the expansion coefficients of the correlation of two circle 

O polynomials, (c) the Fourier coefficients occurring in the cosine representa- 

tion of the circle polynomials, (d) the transient response of a baffled-piston 
acoustical radiator due to a non-uniform velocity profile on the piston. 



O 



X 



^A catalogue record is available from the Eindhoven University of Technology Library 
ISBN: 978-90-386-2290-3 



1 Introduction 

Zernike circle polynomials are extensively used in the characterization 
of circular optical imaging systems with non-uniform pupil functions [1]- 
[8], and, more recently, in the computation of acoustical quantities arising 
from harmonically excited baffled-piston radiators with non-uniform velocity 
profiles [n]^[II]- The circle polynomials were introduced by Zernike [[12] in 
connection with his phase contrast method and the knife-edge test. Further- 
more, they played a fundamental role in Nijboer's thesis [13] on the diffraction 
theory of aberrations, where they were investigated in detail. 

The circle polynomials are given for integer m, n with n — \m\ even and 
non-negative by 

Z-(p,^) = i?H(p)e*'"\ 0<p<l, 0<^<2vr, (1) 

where the radial polynomials Rn are given by 

RtKp) = P""' Pf ""'\2p^ - 1) = E (-1)^ C" ;; ' ) ( ! ) P""-'' ' (2) 



p 

s=0 



,(",/3) 



with Pj^ the general Jacobi polynomial as in [H], Ch. 22 and p = Min- 
imi). It is customary to refer to n as the degree and to m as the azimuthal 
order of Z^. The circle polynomials form a complete orthogonal system of 
functions on the disk < p < 1, with Z™(1,'(9) = e*"*'^, that is, 

i?H (1) = 1, (3) 

and the orthogonality property reads explicitly 

1 27r 

n + 1 



with S Kronecker's delta. In the sequel it will be convenient to set Rn = 
Z™ = for integer values of m, n such that n — \m\ is odd or negative. 

A crucial property of the circle polynomials for diffraction theory is that 
their Fourier transform has the particular simple form 



f f e2"^^^+2"^^^Z^(p,^)di^cip 



1 2n 



27rr 





where we have written u + ifi = pe^"^ and x + iy = re*"^. Equivalently, in 
terms of Hankel transforms (of order m), we have 

1 

J rI-Kp) JM-rp)pdp = (-1)^ ^-^^^^ . (6) 



This formula was given in [12], Eq. (23) and [13], Eq. (2.20). 

By Fourier inversion in Q, using Hankel transforms of order m, it is seen 
that 

oo 

i?H(p) = (-l)'^ j Jn+i{u)J\^\{pu)du, 0<p<l. (7) 



This result, often attributed to Noll [6j, is shown in [12] to follow from the 
discontinuous Weber-Schafheitlin integral, see [13], 15.4.6 on p. 561. The 
integral on the right-hand side of Eq. ([T]) converges uniformly in any closed 
set of p > not containing 1, see Appendix A, and its value for p > 1 is 0. 
Thus the equality in (JTh holds pointwise and not just in an L^-sense. Also 
see [T^, Appendix, Sec. A. 1.1 for a discussion of the result in Eq. ([T]). 

The results in Eqs. (|6|-([7]) are basic to the proof of a number of results of 
the circle polynomials and their radial parts. In [6] , the result in Eq. ([T]) was 
used to derive expressions for the derivative of K^ in terms of i?-polynomials 
of azimuthal orders ?ti ± 1 by employing recurrence relations for Bessel func- 
tions and their derivatives. In [IS]-[in], the two results in Eqs. (|6|-([7]) were 
combined with recursion properties of the Bessel functions to produce the 
scaling formula 

K'iep) = Yl (^n'(^) - K'^'(^)) Kip) = 

n 

n + l 



;E 



n 



n' + 1 



iKt,\{e)-Kt\{e))K{p) 



Here m = 0,1, ... , n' = m,m + 2, ... , and the summation is over n = m,m + 
2,. . . ,n' in which we recall the convention that -R^/"*"^ = Rn'-i — ^ ^^^ ^^^ 
last term in either series. Although Eq. (|8| is normally used for e, p G [0, 1], 
it should be emphasized that they are valid for all complex values of e and 
p by analyticity. The result in Eq. ([s]) is of interest to both the lithographic 
community and the ophthalmo logical community, see [Z], [H], [IS], [16] . 
In [T7], the formula 

n^ip,^) = -^ {l-p'y/'U^{p)e'"'^xlo,i){p) (9) 

n + l 



for the Radon transform of Z^(p,f9) is used to show that 

iV-l 



^rr,/ N J- v^ ^x / /7rfc\ ZTirriK , , 

^n(p) = ^ E^4/^^°^l^)^0^^^' 0<p<l. (10) 



Here m = 0, 1, ... and A^ is any integer > n + m, and Un is the Chebyshev 
polynomial of degree n and second kind. This formula is interesting since it 
gives the R^{p) ioi m < n < N in the form of a discrete cosine transform, 
also see [H]. The result in Eq. ^ was discovered by Cormack in [19], but 
a proof can also be based on the result in Eqs. (|6])-([7]) and the connection 
between Bessel functions and Chebyshev polynomials through the Fourier 
transform, see p3J, 11.4.24-25 on p. 486. 

In the present paper a number of new applications of the results in 
Eqs. ©-(IT]) are presented. A common feature of the problems we con- 
sider is that they all give rise to infinite integrals involving the product of 
three Bessel functions. In Section 2 we consider the problem of finding the 
Zernike expansion of scaled-and-shifted circle polynomials. That is, given 
a > 0, 6 > with a + 6 < 1, we give explicit expressions, involving Jacobi 
polynomials, for the coefficients -ft'^™ (a, b) in the expansion 



Z!^{a + bp'e''>') = J2 ^nT'(a,&) ^„"?'(p'e^''') , < p' < 1 , < # < 27r 



(11) 
see Fig. 1 for the relation between the radial and angular variables of the 
full, centrahzed pupil and the reduced, shifted pupil. 

This result generalizes the scaling formulas in Eq. dsj) which is the case 



with a = in Eq. (11). Furthermore, it gives the analytic solution of the 
transformation problem for the aberration coefficients of an eye pupil when 
the pupil is scaled and displaced. This problem has a long history in the oph- 
thalmological community, see [7] and [20jj for recent work and survey material, 
but no closed-form solution seems to have been found thus far. Moreover, by 
expressing the scaled-and-shifted polynomials as linear combinations of the 
orthogonal terms Z™ one has a handle, via the transformation matrix ele- 
ments K!^^ (a, b), to tackle the important problem of assessing the condition 
of a finite set of circle polynomials when they are restricted to subdisks of 
the unit disk. 

In Section [3] we consider the problem of computing the Zernike expan- 
sion coefficients of the correlation of two circle polynomials. Having these 
expansion coefficients available is of great interest when calculating transfer 
functions in optical imaging. The complex pupil function of the lens is ex- 
panded as a series involving the circle polynomials. The resulting series with 




Figure 1: Centralized, full pupil pe , 0<p<l, 0<'&< 2%, and shifted 



and reduced pupil a + b p' e , < p' < I, < 13' < 27r, where a > 0, 6 > 0, 
a + b<l. 



the appropriate Zernike coefficients represent the amplitude and phase of the 
complex pupil function, including wavelength shift and defocusing. Through 
the analytic expression for the expansion coefficients of the correlation of two 
circle polynomials, one has direct access to the modulation transfer function 
over the full bandwidth of the imaging system. This mathematical device is 
very attractive when calculating, for instance, the pattern transfer in high- 
resolution optical lithography. Hence, when m, n ,m', n' are integers such 
that n — |r7i| and n' — |r7i'| are even and non-negative, and when we denote 
for two functions Z', Z' G L^(]R^) the correlation 

{Z * *covr Z'){U, jJ,) = // Z(z/ + Z/i,/i + /ii)Z*(z/i,/ii)(il/irf/Xi , (12) 



we are interested in finding the numbers T^^J^ such that 



ram'm" lym" ( 1 
n" \2 



4^ ^ nn'n" ^< 



(13) 

Note that Z^ **corr ^™ is supported by the disk around of radius 2 and this 



requires Z^, (gP, "&) at the right-hand side of Eq. ( 13 ). This problem has been 



considered in the optical context by Kintner and Sillitto in [2T], [22] in the 



interest of computing the optical transfer function (OTF) from the Zernike 
expansion of the pupil function. In |2T] , [22] a number of results is obtained 
for the quantities F, but no closed-form solution is given as we obtain here. 



Moreover, the values of the left-hand side of Eq. ( 13 ) are expressed here as 
integrals involving the product of three Bessel functions. 

In optical simulations and/or experiments, a judicious choice of sampling 
points is important. It now turns out that replacing the radial variable p, 
< p < 1, by cosx, < X < 7r/2, has the effect that the variation of the 
radial polynomials Rn is spread out uniformly over the x-range (see Fig. 6 
in Section [T]). This suggests an adequate pupil sampling strategy. In [23] . 
Subsec. 2.4, a matching procedure, using a separable set of sampling points 
(pky'^i) on the disk, for estimating the Zernike expansion coefficients of a 
pupil from its values at the sampling points is proposed. It turns out that 
choosing the radial sampling points as cos x, with uniformly spaced x between 
and 71/2, produces near-optimal results, in the sense that the resulting 
method competes with Gaussian quadrature for all relevant azimuthal orders 
m simultaneously. The fact that the variation of Rn is spread uniformly over 
the x-range suggests, furthermore, to apply this substitution when integrals, 
involving the product of a radial polynomial and a function obtained from 
the pupil function after azimuthal integration, have to be computed. All this 
motivates consideration in Section|4]of the third problem: Finding the Fourier 
coefficients of the radial polynomials Rn in their cosine-representation. Here 
we aim at finding the Fourier coefficients a^^ in the representation 

Ln/2J 

i^lTl (cos x)=J2 <n^2, cos(n -2j)x. (14) 

i=o 

The a™^ will be found explicitly, and from this result it is seen that they are 
all non-negative. It then follows from i?„ (1) = 1, see Eq. (3|, that 

|i?l'"l(p)|<l, 0<p<l. (15) 

This result was proved by Szego, see [21], 7.2.1 on p. 164 and the references 
given there, and the non-negativity of the a™;, was also established by Koorn- 
winder [2S], Corollary 6.2 on p. 113, as was communicated to the author by 
Erik Koelink ^26j . The explicit form of the a^j^ does not seem to have been 
noted before. With the explicit result for the aj^^ available, the computation 
of the above mentioned integrals can be done directly and very explicitly by 
choosing the appropriate sampling points, of the form cos x with equidistant 
X, and using, for instance, DCT-techniques. The fact that the aj^;., being all 
non-negative with sum over k equal to 1, are all small, renders this approach 
intrinsic stability. 
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The use of the Zernike circle polynomials to handle the diffraction inte- 
gral for the acoustical pressure and associated quantities was proposed only 
recently, see P]-[II]. For the case of a circular, flat radiator in an infinite 
baffle, the analytic results of Greenspan |27] on on-axis pressure, far-field, 
power and directivity are generalized and systemized in [9]- [10] while [H] 
deals with radiation from a flexible spherical cap on a rigid sphere. The 
results of Greenspan in [26j, Section VI, on transient responses have been 
briefly considered in the context of the circle polynomials in [2^. Follow- 
ing Greenspan in f2E\, Secton VI, it shall be argued in Section |5] that the 
transient response ^s{t', r) at time t > and at the field point r due to an 
instantaneous volume displacement A at t = of the flexible piston with 
velocity profile v{a), < a < a (=radius piston), is given by 

oo 

<!>s{t; r) = ""^^^fj '^ I Mub)Muw)V{u)udu . (16) 





In Eq. (16), c is the speed of sound, H is the Heaviside function with unit 



step at 0, Vs is the average velocity 



Vs = ^ I I vdS (17) 



7ra^ 



s 



with S the circular piston of radius a, and r is given in cylindrical coordinates 
as 

r = {x, y, z) = (w cos ip, w sin ip, z) , (18) 

while 

b = b{t; z) = VcH"^ - ^2 > (19) 

to be considered for ct > z. Finally, V{w) is the Hankel transform of order 
of f , given by 



V{u)= / Jo{ua)v{a)ada . (20) 



We restrict here to radially symmetric velocity profiles, although many of 
our results can be generalized to the case of non-radially symmetric profiles. 
When now v{cr) is expanded as 

oo 

v{(r) = VsY,^n Rlni^/a) , 0<a<a, (21) 

n=0 



(Zernike expansion for radially symmetric functions on the disk using circle 
polynomials of azimuthal order 0), the result in Eq. ^ shows that 

oo 

^ . r) = g^-^(g^ - ^) J2 (-1)" Un [ Mub) Jo{uw) J2n+i{ua) du . (22) 

Hence, again integrals of the product of three Bessel functions arise. 

A word about the notation. We identify complex numbers z with their 
polar representation p e^^ or their Cartesian representation v + i/i, whatever 
is most convenient in a particular setting. Thus we write things like 

Z^la + p'foe^'''), Z-(z/,/.) (23) 

to denote ^;r(P,^) = I^\p) e'""' in which 

a + p'he'"^' = v + ili = z = pe'^ (24) 

with a, 6 > 0, a + 6 < 1; < p' < 1, ^' G [0, 27r); z/, /i e M, 1/2 + ^2 < ^. 
< p < 1, < ^ < 27r. 

2 Scaled-and-shifted Zernike circle polynomi- 
als 

We shall prove the following result. 

Theorem 2.1. Let a > 0, 6 > with a + h < 1, and let n, m be inte- 
gers with n —\m\ even and non-negative. Then 

Z-(a + 6p'e^^') = Y. KZr'ia,b) ^^^'(p'e*^^') , < p' < 1 , < ^' < 27r , 

n',m' 

(25) 
where for n = \m\, \m\ + 2, ... , n' = \m'\, \m'\ + 2, ... 

K:f{a, b) = TZ^'ia, b) - T^^^^l'^M. b) • (26) 



Here 



rpn 



(-1)^ ^ / Jm-m"{au) Jn"{bu) Jn+l{u)du 




r {<i+p"y- jp-fV- 

{q-q"y. {p+q"y. 



^m-ml'-ua' p(m-m",n"), 

^v-v" ' 



when n — n" > m — m" > 






1—2A'^) p('""~'"'"")/'o d2 



(252-1) 



when n — n" > m" — m > 



otherwise . 



In Eq. (27) we have written 



P 



n — m 



^11 ,^„ii 



^11 I r.^11 



Q 



n + m „ n — m „ n + m 



, P 



Q 



(27) 



(28) 



2 ' 2 ' ' 2 ' ' 2 

Furthermore, P^^' (x) is the general Jacobi polynomial as in [H], Ch. 22 of 
degree k = 0,1, ... corresponding to the weight function (1 — x)'^{l + xY on 
the interval [—1, 1]. Finally, 



1-2^2 



25^-1 



(1 - (a + 6)2)(1 - (a - bY) '^' - (a + 6)(a - 6) , (29) 



1/2 



(l-(a + 6)2)(l-(a-6)2) 
Alternatively, we have 

A = sin a , i? = sin /3 

with a > 0, /3 > such that a + (3 < 7r/2 and 

a = sin a cos (3 , b = cos a sin /3 . 



{a + b){a-b). (30) 



(31) 



(32) 



That is, A and -B can be obtained from the geometrical picture in Fig. 2 
where 7 G [f , vr] is such that sin 7 = a + b. 




Figure 2: Geometric definition of angles a, /3, 7 and side lengtlis A, B from 



a>0, 6>0, a + 6<lin accordance witli Eqs. (31)-(32) and tlie rule of 

sing sin/3 sin 7 -1 



smes 



B a+b 



Proof. By completeness and orthogonality of the circle polynomials, see 
Eq. Q, we have that 



^nT'(a,&) 



n' + l 



1 2-K 



TX 




Z';^{a + hp'e'^'){Z';^\p'e'^')y p'dp'M' . (33) 





We write pe^"^ = a + p'be^^ in which p and ^ are depending on p', d' with 
< p' < 1, < ^' < 27r. Then we get 

1 27r 



^nT'K&) 



ra' + l 



TT 




i?|rl (p(p', d')) e^-''(/'''''') i?J;r'' (pO e-'"'''' p' rfp' rf^' 







We now use Eq. (M) to rewrite Rn {p{p', i!)')) in integral form and change the 
order of integration; this is allowed on account of Appendix A where we show 
that the integral on the right-hand side of Eq. ([T]) converges boundedly for 
all p > and uniformly, to (— 1)^2~ Rn (p), on any set [0, l — e] with e > 0. 
Therefore, 



Kn^ {a,b) = -—- (-1)^ / J„Mu) 



TV 



\ ( 2-K 



J|„l {up{(^, &)) e™'^('''''^') e-'™'''' d^' } ■ R^::'^ (pO p' dp' 



.0 Ko 



du 



(35) 
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Figure 3: Geometric relation between angles x, 7 and side lengths U, V, W 



for Eq. (36) 



We use the addition theorem of Graf, see [m], 9.1.79 on p. 363, 

oo 

Ci{W)e''^= J2 Ci+k{U)Jk{V)e'''r (36) 



fc=— oo 



for integer / and C = J, where W, U, V, x and 7 are related as in the picture 



in Fig. 3. With the variables p', d' , p, i) as in the integral in Eq. (35) and 



displayed in Fig. 1, we use Eq. (36) with 



W = up{p',^'), U = ua, V = up'b, x = ^{p',^'), j = n-^', (37) 
where we note by the comments in [H] after 9.1.79 and 9.1.80 on p. 363 that 



Eq. (36) can be used without any further restriction on U, V, W. Then we 
get 

00 
Jm{up{p',^'))e''-'>^^'-'''^= J2 Jm+k{ua)Uup'h)e'^^--'''K (38) 



fc=— 00 



Hence, in the case that ?7i > 0, 



27r 



MMp', ^0) e'"''(^''''') e-'"''''' dd' 



27r(-l)™' Jm-m'{ua) J-rn'iupb) = 2n Jm-m' (ua) Jm'iupb) 



(39) 
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and so 



oo 

r 

Kfia^b) = 2{n' + 1)(-1)'V^ / Jn+iiu) Jn.-^^{ua) 



/ Jm'{upb)R^\p')dp' 



Lo 



du. (40) 



Using Eq. (38) for m < and noting that J\m\iz) = (— 1)*^ Jmiz) while 



n — \m 



(-i)-^(-ir = (-i)^ 



(41) 



for ?T2 < 0, it is seen that Eq. (40) holds for all integer tu 



Next we use Eq. rt6| with m', n', uh instead of m, ra, 27rr to rewrite the 
integral in [ ] in Eq. ( 40 ) , and we obtain 



^nT'K&)=2(n'+l)(-l) 



n-\-n —m — m 
2 



Jn+lju) Jm-m'{ua) Jn'+l{ub) 

uh 



Using [T^, first item in 9.1.27 on p. 361 

Jn'+l{z) 1 



2(n' + 1) 



{Jn'{z) + j;„/+2(2;)) 



du . 
(42) 

(43) 



it then follows that 



^nT'K&) 



2 



Jn+1 (u) Jm-m' {uo) J„/ {uh) du + 



+ / Jn+l{u) Jm~m'{ua) Jn'+2{ub) du 




,(44) 



and this establishes the equality in Eq. (26) with T's given in integral form 



by the first identity in Eq. (27). The second identity in Eq. (27) follows from 



an application of a result of Bailey, the administrative details of which are 
deferred to Section [H 
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We now list some special cases of Theorem 2.1. 
a = 0. This gives the result of the scaling theory as developed in [T5] . 



also see Eq. (pj). To see this, note that in Eqs. (31)-(32) we have 
A = a = a = , B = b = sin /3 , 



(45) 



and in Eq. (27) only the cases with m = m' give non-zero results. Further- 
more, 



pr'ii) 



,(0,n') 
k 

and, see Eq. (|28l). 



,(0,n') 



b-'pr\2b^-i) = R-^,,ib) 



p + q' = p' + q , p — p' = q — q' 



(46) 
(47) 



since m = m'. Plugging all this in into Eq. (27) yields Eq. (pi). 



b = 0. Only n' = gives non-zero results in Eq. (|27) and then also 

(48) 



m' = (since \m'\ < n'). Now we have in Eqs. (31)-(32) 
A = a = sin a , 



B 



/3 = 



and, see Eq. (28) 



q-V p> = q — 4 1 p — p' = p + (^ 



(49) 



since p' = q' = m' = n' 



0. Thus, when ?Ti > 0, the first case in Eq. (27) 



yields 



^n^pK0)(i _ 2a2)i-(™.o)(-l) = a"^P,(°'™)(2a2 - 1) = RZ+2M = Ki^) 



where we have used that 



P^^'\-^) 



k d('5.7), 



-irp, 



x] 



p(0,m)Q\ ^ ^ 



Therefore, we have in Eq. (25) with 6 = the trivial representation 

Z™(a) = Kia) Z',{p') , < p' < 1 , < ^' < 27r , 
with a similar result in the case that m < 0. 



(50) 
(51) 

(52) 



— a + b = 1. We have in Eq. (31)-(32) in this limit case 

a + /3 = 7r/2, A = sma, B = cosa, 1 -2A^ = 2B^ - 1 = cos2a . (53) 
As a consequence, all T's are non-negative. 



13 



A further interesting observation is that 

ir;^J^'(a, 6) ^ ^ \m'\ <n' <n-\m-m'\ . (54) 

Hence, a Zernike circle polynomials Z"^ can be identified from the set of 
integer pairs (m', n') corresponding to non-zero coefficients when Z^ is scaled 
and shifted. In Subsection 7.1, a detailed computation based on Theorem 2.1 
of the Zernike expansion of the scaled-and-shifted circle polynomials Z^ and 
Zl is presented. 

The validity of Theorem 2.1 will now be shown to extend to all complex 
values of a and h. First assume that a>0, 6>0, a + 6<l. We have from 
Eqs. Q-Q that 

^r(pe^') = {pe^'')\-\ pS^(2p' " 1) > (55) 

2 

where ± = sgn(m). We have from 

pe^'"" = a + bp'e^'''' , (56) 

p2 = a^ + b\p'f + 2abp' cos ^' (57) 



that 



n+|Tn| 



Z„"^(a + bp' e'^) = Yl P^(P 5 «' ^) e^^"" ' (58) 

where Pk{p' ', d, b) depends polynomially on p', a and b. On the other hand. 



from Eqs. (29)-(30), we have that for any polynomial p 

p(l - 2A^) p{2B^ - 1) = p{-x + y)p{-x - y) (59) 

is an even function oi y = ((1 — (a + &)^)(1 — (a — 6)^))-*^/^ for any value of 



X = {a + b){a — b). Consequently, the right-hand side of Eq. (59) contains 



only even powers of y. We conclude that any of the T's considered in Eq. (27) 



depends polynomially on a and b. Hence, the relation in Eq. (25) extends to 
all a, 6 G C by analyticity. 

An important consequence of this extension is that now also the trans- 
formation matrices {K^^ {—a/b, 1/b)) corresponding to the inverse transfor- 
mation z H> —a/b + z/b can be considered. Accordingly, when the degrees 
n, n' are restricted to a finite set {0, ...,A^}, the matrices corresponding to 
z ^^ a + bz and z H- —a/b + z/b are each other's inverse. Having expanded 
the shift-and-scaled circle polynomials in terms of the orthogonal functions 
Z^ , we have now the opportunity to deal with the problem of assessing the 
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condition of the set of circle polynomials of maximal degree A^ as a linear 
system when they are restricted to an arbitrary disk in the plane. Indeed, 
the condition number is given as the square-root of the ratio of the largest 
and smallest eigenvalue of the Grammian matrix, and this Grammian ma- 
trix and its inverse are expressible in terms of the appropriate transformation 
matrices {K^^ {a,b)). Such an effort is already worthwhile for the case of 
restriction of circle polynomials to a disk p < e with e < 1 (pure scaling), 
and the author has found for this case useful and simple estimates for the 
magnitude of these condition numbers. This case is much simpler than the 
general case since the transformation matrices decouple per m. 

3 Zernike expansion of the optical transfer 
function 

In this section we assume that we have expanded a pupil function P{p,^) 
(vanishing outside p < 1) as 

nP,^) = E^^^^^^^('^^*'')' 0<P<1' 0<^<27r, (60) 

z — / 7j- 
n,m 



m 
n ■ 



with coefficients 7™ that can be obtained by using orthogonality of the Z, 
Writing u + ifi = pe^^ with z/, /i G M and identifying P{v^p) = P{p^d), 
compare end of Section [T| it is required to find the Zernike expansion of the 
OTF (optical transfer function) 

(P **corr^)('^,/i) = // P{u + h'i,p + pi)P*{yi,pi)dvidpi (61) 

that vanishes outside the set z/^ + p^ < 4, see Fig. 4. Thus, considering the 



expansion in Eq. (60), it is required to compute (the Zernike expansion of) 
Z^ * *corr Z^ for integer n, -m, n\ m' with n — \m\ and n' — |r7i'| even and 
non-negative. We maintain the p, g-notation of Eq. (28). 

Theorem 3.1. We have, with u + ip = p e^^ and < p < 2, 

(Z-**eo.r^„"?')(z^,/i) = 

oo 

J u 
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and 



(n + l)(n' + l) 



m' \/ i'd\ 



TT^ 



;zr*w^„"?)(pe 



V^ n + i -f^mrn'm" rym" / 1 ^ J-d-\ 



A-K 



-pmrn'm" ym" /"!„„* 
-■- nn'n" ^n" Kof^ 



where F is non- vanishing for m" = m — m' only, and in that case 



(63) 



-pmmm 



]{n + l){n' + l){-iy 



Jn+liu) Jn'+l{u) Jn"+l{2u) 



du 



u^ 



(64) 




Figure 4: Schematic representation of the autocorrelation function P * *corr P 
of a non-uniform pupil function P as an integral over the common region of 
two disks of unit radius with centers at and pe*'', respectively. 
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Proof. We have by Parseval's theorem and Eq. ([s]) that 



(zr*w^„"?')(z^,/i) 



7[Zl\u + ■,^l + ■)]{x,y){:F Z^'y{x,y)dxdy 



II e~'-"''''-'^y{J'Z::^){x,y) {:FZ:^'y{x,y)dxdy 



oo 27r 




Jn+i(27rr) J„/+i(27rr) _ 



g --^ipr cos{^-ip) ^i(m-m')ip ^ ^^ i 







(65) 



where we have written x + iy = r e^'^ and where we use that ux + fiy 
p r cos{'& — ip) . Now 



27r 



'2TTipr cos{'d—ip) ^ikip 



e'"^ dip = 2TT{-i)'' e'^"" Jfe(27rpr) 



(66) 



and inserting this into Eq. (65 ) with k = m—m', noting that i" " (— i) 
(—1)^"^ and substituting u = 27rr, we obtain Eq. (62). 



see Eq. ^, that 



Next to show Eq. (64), we note from orthogonahty of the Z™, {^pe^^^ 



"n" V2^ ^' 



ri7„r = (^ + i)y + i) / / (z™ **,„,, z™')(z/,/.)(z™"(iz.i/i))*rfz.d/i 



TT" 



(67) 



Again using Parseval's theorem, together with 
T{Z **,,,, Z')=TZ-{TZ'r , {J' Z"{\v,\p)){x,y) = A{J'Z"){2x,2y) , 

we obtain 



(6^ 



y'mmm 



A{n + l){n' + 1) 



TT^ 



(^Z™)(x,y)(J-Z;:?')*(a;,i/)(^Z^;')*(2x,22/)dxd2/. 

(69) 
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Then inserting Eq. (pi) and using polar coordinates x + % = r e*"^, we obtain 

oo 

r-^y = ^-t±^M±}l,n-n'-n" ^ ^^^^(2vrr) J„,,,(2vrr) J„„,,(4.r) ^ ■ 



27r 



The proof is completed by the substitution u = 27rr. Also note that n — n' 
and n" have the same parity when m" = m — m' . 

Notes. 

1. We have by Eq. Q that 

(^r * *co„ ^„"j')(^, /^) = ^ (-1)^-^' e^^"'^'' [Q--™; + Q:r+Tn' +2] (71) 

and, for the case that m" = m — -m', 
where 

oo 

Q^j{a,b,c)= / Ji{au) Jj{hu) Jk{cu)du , (73) 



and where the Q's in Eq. (71) are evaluated at (a = 1, 6 = 1, c = p) and the 



Q's in Eq. (72) are evaluated at (a = 1, 6 = 2, c = 2). 



2. We shall use in Section p^ a result of Bailey f28| to show the follow- 
ing. We have Q"„^^ = when n" < n + n\ and when n" > n + n', we 
have 

^ (IK + n + n')y- (IK - ^ - ^0)! A.n+n'+l p{n,n') ..^ p(n',n) . ^ 

(1(^// + ^ _ n/))! (i(^// + n' - n))! 4^ ^B^^^^iuj z^^^^^^^iuj , 

(74) 

also see Note 2 at the end of Section 4. 

3. From a result of Bailey |28| it also follows that Q"^, and Q"„',+2 
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vanish when p > 2, but this is aheady clear from the fact that the Z's are 



supported by the unit disk. The result in Eq. ( 71 ) takes a more complicated 



form when < p < 2, see Section |6] for more details. 

4. In Subsection 7.2, a detailed computation, based on Theorem 3.1, for 
the OTF corresponding to P = Zq is presented. 

4 Cosine representation of the radial polyno- 
mials 

We shall prove the following result. 

Theorem 4.1. Let m, n be integers > with n — m even and non-negative. 
Then 

Ln/2J 

i?™(cosx) = y^ a„_2j cos(n — 2j) X , (75) 

j=0 

where for integer A; > with n — k even and non-negative 

ak=e,^^_{\)\Pl^''\m ■ (76) 

Here Sq = 1, ei = €2 = ... = 1 (Neumann's symbol), and 

n—l n+l n—r n+r l—r l+r 

^=^-' ^=^-' ^=^-' ^ = ^-' ^=^-' ^=^-' 

(77) 

where 

/ = max(?7i, /c) , r = min(?7i, k) . (78) 

Proof. We have from Eq. (ul) 



00 



K!^{cosx) = (—1) 2 / J^_^_i(u) Jm{u cos x) du . (79) 


Next we note that 

R';^{cosx) = icosxrP^^\cos2x) (80) 

2 

has non-vanishing Fourier components bk e^^^ only for integer k of the same 
parity as m. For such k we shall show that 

27r 
/ JmjU COS X) e^^^ dx = J7n-k(^u) Jm + k(^u) . (81) 

2tt J 2^2^ 
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Indeed, abbreviating "the coefficient of e*™^ in" by Cm, we have by the gen- 
erating function 



e''''''y= J2 M^)e 



imy 



12) 



that 



27r 
1 /" 1 

-— / Jm{u COS x) e*''^ dx = --C„ 

ZTT J ZTT 





27r 



iu cos X sin y ikx j 



2-K 



Or) 



27r 



,|msin(a;+y) g-|msm(x-i/) gifcx ^^ 



27r 



Gn 



Gt, 



27r 



5^ J^A\u)JmA\u) j e*'"^ 



(x+y)-im2{x-y)+ikx i 



m,i,m2=— oo 



/ ^ ^mi I, 2^/ '^"2 V 2^' ^ 



+7712)?/ 



mi , 7?^2 — — oo, 
m2 — rrii — fc 



Jm-k{\u) Jrn±k{\u) . 
2 2 



^3) 



Then we have at once from (79) that 

2-K 



K^{cosx)e'-^'' dx = (-1)"2™ / J„+i(m) J„^-k(lu) Jm+k(lu) du . 



1 

2^ 



When A; > and ?7i — /c > 0, we have then 



27r 



1 

2^ 



R^{co^x)e^^^dx = 2{-lY-^Q:^^{lX2) 



lA) 



Ih) 



see Eqs. (73) and (74). For this case, the result follows from Eq. (74), the 



fact that Pf'^'\Q) = {-1)P Pi^'^\0) and the fact that a^ = a_k- When 
A; > and m — k < 0, the result follows in a similar manner, by using that 



Jm-k(^u) = (— 1) 2 J k-m (^u) and Eq. (74) together with some easy admin- 
istration with signs. 
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Notes. 

1. It is straightforward to generalize the resuh of Theorem 4.1 to the repre- 
sentation of R^iycosx) with < f < 1. Now Q^'ti „+fc (f , f , 2) appears in 



Eq. (85) and a^ in Eq. (76) becomes 



2 ' 2 



ak{v) = eu^-^{\v)' Pi,'<''\x)Pi,'<^'\-x\ (86) 

where a; = (1 — 1>^)^/^. Also, see Section [6J and Eq. (27). 
2. In using this result it is convenient to note that 



see [S], 22.3.1 on p. 775. 



3. In Subsection 7.3, a table is presented in which the radial polynomi- 
als i?^, integer n, m with < n, m < 8 and n — ?7i > and even, are given 
in the polynomial representation of Eq. (|2]) and in the cosine representation 
of Theorem 4.1. 

5 Transient response from a flexible, baffled, 
planar, circular piston 

In this section we consider sound radiation from a circular, baffled-piston 
radiator with (possibly non-uniform) velocity profile f (r^) that vanishes for 
Ir^l > a, where r^ is a point (x^, |/s, 0) in the plane z = Q and a is the radius 
of the piston. The complex amplitude p(r, u) of the radiated pressure at a 
field point r = {x,y,z) with z >0 due to a harmonic excitation exp{iut) is 
given by Rayleigh's integral or King's integral (in the case that v is radially 
symmetric) |29], [30] 



? On Cn i P ^ 

P(r,c^) = -^^ / v{Y,)—^dS 



s 



°° ^(«2-fc2)l/2 



ipock / -^———j^Jo{wu)V{u)udu. 
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Here po is the density of the medium, c is the speed of sound in the medium, 
k = u/c is the wave number with u the radian frequency, 5* is the disk of 
radius a centered at the origin and perpendicular to the 2;-axis, r^ is a point 



{xs, Us, 0) = {a cos f?, a sin t9, 0) 
with < cr < a on S", r is a field point 

(x, y, z) = {w cos V', wsmip, z) 



59) 



(90) 



with ti; > 0, z > 0, and r' = |r — r^l is the distance between r and r^, see 
Fig. 5 for geometry and notations. Furthermore, {v? — /c^)^/^ is defined to be 
i y/k^^-v? for < u < A; and \/v? — k^ ioi u > k with non-negative ^, and 
V{u) is the Hankel transform of order 0, 



V{u) = / Jo{ua)v{a) a da , u>0 





(91) 



of the radially symmetric profile f (Ir^l). 

Bouwkamp has generalized the King integral representation to the case 
of non-radially symmetric profiles 

oo 

v{a,^)= J2 v"'{(^)e'"'^ , 0<(T<a, 0<^<27r, (92) 

m=— oo 

as follows. There holds El] 



p{r,uj) =ipock ^ 



,-2(«2-fc2)l/2 



m=— oo 



,, _ p.i/2 -Jmiwu) V^'iu) u du e*"^'^ , (93) 



where V"^ is the Hankel transform of order m of v"^' 



V^iu) 



Jm{ua) v"^{a) a da . 



(94) 



In defining the transient response ^sit; r) at time t > and field point 
r due to an instantaneous volume displacement A at t = 0, we follow 
Greenspan [27], also see Harris [32]. Thus, the transient response is given by 



$5(t;r) 



vra^ K 



£-1 



_ipoUJ 



it) 



(95) 
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r = {x,y,z) 




Figure 5: Geometry and notations for baffled-piston sound radiation, where 
a is the radius of the piston and 



^s = ixs,ys,0) = ((T cos ^,cr sin ^,0) , 

r = {x,y,z) = {w cos ip,w sin ip, z) , 

r' =\r- r,| = {{x - x,)^ + {y - ysf + z^f'^ 



where Vg denotes the average speed 



(96) 



and £ is the inverse Laplace transform that is apphed in Eq. (95) to the 



function p(r, ijj)/i po ^ in which the wave number k = u/c is replaced by s/ic 
with s the Laplace variable. Greenspan in [27], Section VI uses the King 
integral (for the case of radially symmetric v), together with an identity for 
the inverse Laplace transform of m~^ exp{—zm), m = {v? + s^ jc^Y^'^-, to 
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obtain 



cAHict — z] 
^sit ; r) = \— — ^ / Jo{u a{t ■ z)) Joiuw) V{u) u du , (97) 





where 



a = a{t- z) = VcH^ - z^ , ct>z, (98) 

and H{x) = 0, 1/2 and 1 according as x < 0, x = and x > (Heaviside 
function). Using Bouwkamp's generalization of the King integral represen- 



Eq. (95) 



tation of the pressure, see Eq. (93), this yields for the transient response in 



cAH{ct- z) 



Y^ Jo{ua{t; z))J,r, (uw) V'^{u)udu- e'""^ . (99) 



m=— oo 



We now develop the general velocity profile v{rs) = v{a,'d) into circle 
polynomials according to 

v{a,^) = VsJ2 <RtKa/a)e'^'^ , (100) 

n,m 

SO that v^{(t) is given by 

i;™(a) = K Xl""^n'"'(^/") ' 0<a<a. (101) 

n 

Noting Eq. (|6]), so that 

a 

I Rt\cy/a) Jmiua)ada = (-1)"^ - Jn+iiua) , (102) 

J u 



we thus get that 

Mt; r) = "^^^"'~"^ J2i-iy-^u-Q:,^,^^{a,a,w)e-^ , (103) 



vra 



where 



oo 

Qn+i,m(c^)«)^) = / Mua) J^{uw) Jn+i{ua) du , (104) 
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with Q as in Eq. (73). We shall consider the integrals at the right-hand 



side of Eq. (104) further in Section 6 Until now, except in Eq. (62), the 



parameters a, 6, c occurring at the right-hand side of Eq. (73) were such 



that they do not occur as the sides of a triangle. In that case it follows from 



Bailey's result in Eq. (110) that Q^j^^^^a^a^w) = with integer n and m 



such that n — \m\ is non-negative and even. In (104) we must now allow also 
a,a,w > that do occur as sides of a triangle. 

There are some special results for the case of radially symmetric profiles 
V. Greenspan in [27|, Section VI, has shown that 

A 

<l>s{t;r) = - ^ \~ ^' vUw'^ + a'^ -2wa cos aY^'^] da , (105) 

vr"' a^ K / 



where 







A 



arccosi - 

n 



fw'^ + a"^ — a^ 



2wa 



a < \w — a\ , 

\w — o-\ < a < \w + a\ , 
a > \w + a\ , 



(106) 

holds for velocity profiles v^^\a) that are a multiple of {1 — {a / a)'^)"' H {a — a) , 
with H the Heaviside function as before. In [33] it has been noted that 



Eq. (105) holds for general radially symmetric profiles v. Thus with 

oo 

v{^) = VsYl < ^2n(^/«) , < a < a , 



n=0 



there follows 



n=0 



W] 



Txa 



(107) 



(108) 



where, see Eqs. (103)-(104) 



(-l)"Q2n+l,o(^>«>«^) 



— / RlJ^iw"^ ^-a"^ -2wa COS af^/ a] da 
Tca J 



1 
Txa 



PJ2 



w"^ + a'^ 



1 cos a da 



(109) 
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The last expression in Eq. ( 109 ) is further elaborated in [33] , using the addi- 
tion formula for Legendre polynomials, to yield a finite series expression in 
terms of Legendre functions and sine. 

6 Infinite integrals involving the product of 
three Bessel functions 

We consider in this section the integrals 

oo 

Ax^u{a, b,c) = / Jx{ax) J^{hx) Jy{cx) dx , (110) 



where a,b,c > and A, /i, u are non-negative integers. There is quite some 
literature on these integrals and more general instances of them, see [31], 
Sec. 13.46 on pp. 411-415, [28], [35], Sec. 19.3 on pp. 349-357, [36], Sec. 13.4.5 
on pp. 331-335. These general results can become quite unmanageable; it is 
the point of this section that, in the special cases that we consider, often con- 
cise and manageable results appear, often in terms of the radial polynomials 
themselves. 

For the case that a + b < c (and a > 0, 6 > 0, c > 0), the integral in 



Eq. (110) is given by Bailey, [2H], Eq- (8.1) (with a minor correction in which 
the c^+^+^ in the denominator at the right-hand side should be replaced by 

c), as 

oo 

J\ {cu sin a cos (5) J^{cu cos a sin /3) J^ (cm) du = 

r(|(l + A + /i + v)) sin^ a cos^ /3 cos'' a sin^ /3 
cr(A + 1) r(/i + 1) T{\{1 - A - ^ + I/)) 

■ 2i^i[|(l + A + /i-z/),^(l + A + /x + z/); A + 1; sin^a] ■ 

• 2Fl[\{l + \ + ^i-u),\{l + \ + ^i + u)■, /i + 1; sin2/3] , (111) 

where a, /3 > are such that a + (3 < -n 12. In Eq. ([27J) the choice 



c = 1 , a = sin a cos /3 , h = cos asin/3; \ = m — m' , fi = n' , u = n + 1 

(112) 
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is made. When m — m' > and A = sin a, B = sin/3, we have 



Jm-m'{au) Jn'{bu) J„+i(m) du 



r(|(n + n' + m - m') + 1) a""-""' h""' 
V{m — m' + 1) V{n' + 1) T(^{n — n' — m + m') + 1) 



iF2{i2{n' — n + m — m'), i^{n + n' + m — m!) + 1] m — m' + 1\ A) ■ 



iF2{i^{n — n + m — m), :^{n + n' + m — m) + 1; n + 1; B 



:ii3) 



Recall that both n — m and n' — m' are even. Hence n — n — m' + m' is 
an even integer, and when this even integer is negative the whole expression 
(113) vanishes dues to the r(|(n — n' — m + m') + 1) in the denominator. 



Using [14j, 15.4.6 on p. 561 and 22.4.1 on p. 777, we have for j = 0, 1, ... and 
7,5 >0 



iF2(-j,7 + l + 5 + j;7 + l;a;) 



^•^(^ + ') P^^)(l-2x). (114) 



r(7 + i+j) ^ 

Therefore, 

iF2{\{n — n + m — m), ^{n + n + m — m) + 1 ; m — m + I; A^) = 
= iF2{^{n — n + 771 — 7n'),m — m + 1 + n + 

+ ^{n — n — m + m!) ] m — m! + 1; A^) = 

U(n' — n — m + m'))\ (m — m')\ (m-m' n') , .o. , . 

din - n' + m - m'))\ i(n-n'-m+m')V ^' V ; 

and similarly 



iF2(l(n' — n + m — m), l(n + n' + ?7i — m) + 1 ; n' + 1 ; B^) 



{Ijn-n' -m + m')y.{n'y. ^(„' ,„_„') 

{l{n + n' -m + m'))\ \{n-n' -m+m') 

(7,5) 



'l-2fi^ 



(116) 



Using that Pf'°\—x) = (— 1)-' PJ ' (x) in Eq. (116) and some further ad- 
ministration with F-functions and factorials then yields the first instance in 
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Eq. (27). In the case that m'—m > 0, we use that J„ 



4)r 



Jrr 



and so we can apply the first instance in Eq. (27) with m and m' interchanged. 



This requires a careful administration with q, p, q', p' as well as with the signs 
(— 1)^~P in Eq. (27). Doing so, the second instance in Eq. (27) follows. 



The result in Eq. (27) has been proved now for the case that a + b < 1. 

, a > 0, 6 > 0, follows by taking the limit case 



However, the case that a+b 



in Eq. (27) and observing that the integral in Eq. (27) converges uniformly 
in (a, 6) G [e, 1] x [e,l] for any e > since Jkiz) = 0{z~^^'^), z -^ oo and 
Jk{z) is bounded in z > 0. 



We shall now consider the Q-integrals in Eqs. (73) and Eqs. (84), (85) 



These can be treated in very much the same way as the integrals in Eq. (113) 



that arise from Eq. (Ill) by making the choice as in Eq. ( 
we have here the limit case a = b — ^ 



Q 



n"+l, 



1,1,2) 



IQ 



n"+in 1 



^c. As to Eq. (|74|), we note that 

n + 1) 



12). Note that 



by (n, n', n" + 1) and take 



1). We thus need to replace {m — m' 



A' = B' 



n 



(117) 



in Eq. (113). In particular, the Q-integral vanishes when n' + n — n" < 0, and 



when n' + n — n" > the 1F2 that arise in Eqs. (115), (116) should be written 



down with the replacement just mentioned and the choice in Eq. (117). This 

then yields Eq. (74). 

similar fashion by replacing {m — m',n',n + l) 

with a, b, A, B as in Eq. (113). 



Next, the Q- integral in Eq. (|84|) can be handled in a 

k),\{m + k),n + l) 



by {\{m 



We now consider the integral in Eq. (110) for the case that a,b,c > 



while non of the numbers a, b, c is larger than or equal to the sum of the 



other two; this occurs in Eqs. (62), (71), (73) when < p < 2 and in Eq. (104) 



where no other restriction on a, w, a is imposed than being non-negative. 
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Here the result in [31], Eq. ([T]) on p. 413 can be used. We have 



J\{u cos ip cos ^) J/i(Msin(y9sin$) J^{u cos ^) du 



cos''* if cos'^ $ sin'^ yj sin'^ $ cos'' d 



r(z/ + i)(r(/. + i))2 

2^(-l) (A + /X + 2A. + 1) ^ 



fc=0 



fc!r(A + A; + l) 



r(^(A + /i + t/ + i) + fc) 
r(^(A + /i-i/ + 3) + fc) 



^,z^ — 1 — A — u ,A + u + i/ — 1 , _, _, o„ 

Fi ( - k, ^ + k + l;u + l] cos^ i3 



2 2 

• 2Fi{-k, X + fi + k + I; /i+1; sin^ </)) • 

• 2Fi{-k, X + fx + k + l; fx + l; sin^ $) 
when Re(A + /i + z/) > —1 and cos-;? 7^ ±cos($ ± ip). 



:ii8) 



As to the integral in Eq. (104), we use Eq. (118) with \ = 0, jj, = m 



u = n + 1 with n, m non-negative integers with n — m even and non-negative 
(the case that m < will be considered afterwards). We then have 



/: = 



Jo(m cosy? cos $) Jm('U sin;/) sin $) J„+i('Ucos'i9) du 



sin™- (p sin™- $ cos'^^-^ 1!} 
r(n + 2)r2(m + l) 



00 
^(-l)'=(m + 2fc + l)(^ 



fc=0 



r(m + fc + l)\2 



r(A; + i: 



r(|(n + m) + fc + l) 
r(-^(n-m) + A; + l) 



■ 2i^i ( ^ A;, — — + /c + 1 ; n + 2 ; cos^ ^ 

,2 _\ 17 ^ L _ , L , 1.™ , i.„- 2, 



2-Fi(— A;, m -h A; + 1; m -h 1; sin yj) 2-Fi(— A;, m -|- A; -h 1; m -|- 1; sin $) . 

(119) 
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Next, we use (114) to write the product of the two last 2-^1 in Eq. (119) 



multiphed by sin™ (p sin™ $, in terms of Jacobi polynomials as 

sin™ ^ p^™'0)(l - 2 sin" (^) sin™ $ P^(™'0)(1 - 2 sin^ $) . (120) 

By letting a = | — ^9, /3 = | — $, this latter product can be written as, see 
Eq. (§, 

i?™+2,(cos«)i?™+2,(cos/3), (121) 



where it also has been used that P^™"' [x] 



-1)'' Pf^ '"^(x). Furthermore, 



we note that the terms of the series in Eq. (119) vanish when k < ^{n — m) . 
Letting p = ^{n — m), writing k = p + r with r = 0, 1, ... , and noting that 
m + 2k = n + 2r, we see that 



/ = 5^ (-l)P+''(n + 2r + 1) 



r=0 



r! (n + 1) 



cos"+^ ^ 



■ 2Fi{—r, n + r + 1; n + 2; cos^ f?) ■ 
•i?™+2,(cosa)P™+2,(cos/3). (122) 



We shall next show that 



-lY{n + 2r + 1) 47 7^ a;"+^ 2i^i(-r, n + r + 1 ; n + 2 ; x') 

r! m + 1 ! 



-"^n+2r+l(^) -^■' 



n.+l 



n+2r-l 



,X 



(123) 



where we have -R^li = for the term with r = in the series in Eq. (122) 
As a consequence, it follows that 

^ ~ Joiu sin a sin l3) Jmiu COS a COS l3) Jn+i{u COS ■&) du = 



= i-iy Yl (KlUi(cos^) - K+^,_,(cos^)) i?:r+2.(cosa) i?:r+2.(cos/3) 

r=0 

(124) 
under the condition that cos(/3 ± a) = =1= cos($ =1= y?) 7^ cos -(9. 
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To show Eq. (123), we note that for r = either side of Eq. (123) equals 
^n+i ^y ^^, 15.4.1 on p. 561. For r = 1,2,..., we use 



?n+l 



^n+2r+l^^ 



(x) = x"+^ Pj°'"+i)(2x2 - 1) = (-1)'' x"+^ P,("+i'°)(l - 2x') 



r! [n + Ij! 



and 



pn+l 
-'^n+2r 



r-1 „n+l 



.i(x) = (-ir^x^ 



(125) 



2i^i(-r + l,n + r + l; n + 2; x^) . 



(r-l)!(n + l)! 
Then the identity to be verified becomes 

{n + 2r + 1) 2Fi{-r, n + r + 1; n + 2; x"^) = 

= (ra + r + l)2-Fi(-r,n + r + 2; n + 2; x^) + 

+ r 2Fi{—r + l,?T, + r + l; n + 2; x 

and this is the contiguity relation [H], 15.2.14 on p. 558 with 

b = n + r + 1 , c = n + 2 , 



(126) 



, X j , 



a = —r 



z = X 



(127) 



:i28) 



Hence Eq. ( 123 ) has been estabhshed. 



The result in Eq. (124) has been shown for the case that m, > 0. When 



m < it holds as well when m on the right-hand side of Eq. ( 124 ) is replaced 
by \m\ while p is maintained to be ^(n — m). This follows from Eq. (124), 
with \m\ instead of m, and the fact that Jm{x) = (—1)™ J\m\{x). 
Choosing A = 0, /i = n + 1, z/ = m in Eq. (118) yields 



Jo {u sin asm (3) J„+i {u cos a cos (3) Jmiu cos 'd) du 



-1)" E (<2.(cos^)-Pi"^2.+2(cos^)) K:U(cosa) K:1+2.(cos/3) . 

(129) 



fc=0 



The choice {A, /i} = {m, n+ 1}, z/ = in Eq. ( 118 ) does not give an appealing 
result in terms of the radial polynomials, also see below. 
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An important special case of Eq. (124) occurs when cos -(9 = 1. Since 
i?;^t}(l) = 0, RlXlr+iil) = 1, r = 0, 1, ... , only the term with r = does not 
vanish, and there results 



oo 

/ 



Jo('Usinasin/3) Jm(-ucosacos/3) Jn+iiu) du 



'iyK^{cosa)K^{cosf3) 



(130) 



This can be shown to be equivalent with the result in Eq. (27), upper or 
middle case, with m = m' . Furthermore, the choice sin a = (so that 
Jo(Msinasin/3) = 1) can be shown to yield Eq. ([7|). 



We finally consider two more choices in Eq. (118). The first choice has as 



a special case the result of Eq. (27), and the second choice is of relevance for 



the computation of the Q-integrals that occur in Eq. (|71|) for the evaluation 
of Z"!^ * *corr ^™ at a particular point u^ jji. 



Thus, with A = ra, /i = Z, i^ = n' + l in Eq. (118) where n, /, n' are integers 



such that ra, n' > while n — n' and / have the same parity we get in a similar 



manner as Eq. (124) was proved the following result. There holds 



Jn(Mcos(y9cos$) J;('usin(/)sin$) J„/+i(-ucos'i?) du 



n — n -j-Z \ r) / \ I/I 

^-1) 2 (cos(y9cos$) (sinv9sin<l>)' ' 



fc=max{0,|(n-n'-|/|)} 



n + k + 

\l\ 



>(",! 



y{n,\l\) 



k+\l\ 



Pr''Hl-2sin»Pr'inl-2sin^$) • (131) 



\{n' 



n 



> is 



For the case that ■(9 = 0, only the term with k 

non- vanishing, and this then yields the result of Eq. (27) with I = m — m' 

and n and n' interchanged. 

Next, with X = n, fi' = n' + 1, u = I in Eq. (118) where n, I, n' are integers 



such that n, ra' > while n — n' and / have the same parity and n + n' > 
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we get the following result. There holds 

Jn{u cos ip cos ^) Jn'+i{u simp sin ^) Ji{u cos i!}) du 



'-!)" 2 (cos99cos$)" (sinv3sin$)"'"^^ 



fc=0 



n + n' + k + 1 



n + k 
n 



p,n -,.,n,^^ _ 2gj^2 ^) pyn +x,n;^^ _ 2 gj^2 ^) 



(132) 



This vanishes in the case that cos"!? = 1. For the Q-integrals in Eq. (71) this 



gives a series result by taking in Eq. (132) 



cos if = cos$ = siny9 = sin$ = l/v2 , (133) 

replacing ^u in the integral by u and setting 2 cos d = p. 

7 Examples 

In this section we present worked out examples of our main results. 

7.1 Examples for Section [2] 

We use Theorem 2.1 for the computation of the Zernike expansion of the 
scaled-and-shifted circle polynomials Z^ and Z^. 

7.1.1 Computation for Z^ 

We have 

Z°(a + 6p'e^^') = Y. i^r ^."?'(p'e^''') , (134) 

7i',m' 

where we have K^^ Z™ ^ only if n' and m' have the same parity and 



m'\ <n'<A — m', see Eq. (54). This leaves us with the cases 



\m'\=0, n =0,2,4; \m'\ = l, n' = l,3; \m'\=2, n=2. (135) 
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-0',-m' 



Furthermore, K^^, = K^!^ and so it is sufficient to do the computations 
for the cases m' = 0, 1, 2 in Eq. (135). 



m' = 0. a. K^^ = T°° - T^" by Eq. (g with T™ and T°° given by Eq. (|27|), 
middle case, as 



TOo> = q = 2-p' = q' = 0]= P^'il - 2A') Pfi'^B' - 1 



(0,0), 



i2n p(0,0). 



2 p(0,2). 



l2^ p(0,2). 



(136) 
TO^V = g = 2 ; p' = g' = 1] = 62 pj'^'-^(i _ 2A^) PI''''>{2B^ - 1) . 

(137) 

m' = 0. b. K^^ = T4°2° - r4°4° with T^^ given in Eq. (136) and T4°4° given as 

T^^{p = q = 2-p' = q' = 2] = b' pI,'''\i - 2A^) pI,'''\2B^ - 1) . (138) 



m' = 0. c. K^2 = ^44 - ^46 with T47 given in Eq. (138) and T^^ = 0. 



m' = 1. a. K^l = T^l - T^^ by Eq. (g with T^^ and T^g^ given by Eq. (|27|), 
middle case, as 



(1,1)/ 



\2\ d(1.1). 



T^l[p = q = 2;p' = 0,q' = l] = labP^'''>{l - 2A') P^'''>{2B' - 1 



(139) 



3 p(l>3)/ 



l2\ p(1.3). 



T4"31P = g = 2; p' = l,g' = 2] = Aab' P;,'''>{1 - 2A') P',''-'\2B' - 1) . 



(140) 



ni 



' = 1. b. K^^ = T4O3I - T4O5I with Tfgi given by Eq. (139) and T^^ = 0. 



m' = 2. a. K^l = T^i - T^l by Eq. (|26| with T^^ given by Eq. (|27|), middle 

CclSG, 3iS 

2,2 0(2,2)/^ ^ .2n 0(2,2), 



T4O2 = [p = g = 2 ; p' = 0, g' = 2] = 60^52 p^^'''\i _ 2^^) Po^2-2)(2ij2 _ ^^ 



:i4ii 



and T4I = since n — n" = Q < 2 = \in — 'm"\, see Eq. (27). 



There remains to be calculated the right-hand side of Eqs. (135)-(140) 
with P^^' the Jacobi polynomials and 1 — 2A'^ and 2B^ — 1 given in terms 
of a and b by Eqs. (|29|)-(30). In general, one can use that 



Pt''\^) 



(A; + 7)! 



k\{k + -f + 6)\ 



E(n '':!"Tr" (^-i)' 



1=0 



2'(/ + 7)! 



(142) 
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>(7,5), 



\k p(5,7)/ 



together with Pj: (— x) = (— 1) P^ (x). For the present purposes it is 
sufficient to know that 



p{0,0)| 



x] 



3^2 
2' 



,(0,5) 



x^-l; Pr'^(x) = (l + ^(5)x-i(5; Po^^ Ha;) = 1 • (143) 



Using Eqs. (143) and (29), (30) in Eqs. (135)-(140) yields 



t: 



00 
40 



6a^ + 66^ + 2Aa'b^ - 6a' - 6b 



T^^ = b'{12a' + 4b'-3) 
Til = lab{8a^ + 8b^ - 4) 



t: 



00 

44 

01 . 



T- = 4ab' , 



Til = 6a'6' . 



(144) 
(145) 
(146) 



This then gives 



K' 



k: 



00 
40 

00 
42 



,2t2 



6a^ + 26* + 12a^6^ - 6a' - 3b' + 1 



36* + 12a'6' - 36^ 



K^l = 12a% + 8ab^ - 6ab 



Ml 



K'i = 6a%^ . 



/^OO _ ,4 



K°3^ = 4a6=^ 



(147) 
(148) 
(149) 



Hence 



Z°(a + bp' e'""') = [{6a^ + 26^ + 120^62 - 6a^ - 36^ + 1) Z° + 
+ (36^ + 120^62 - 36^) Z° + b^Zl] + [(120^ + 806^ - 6ab) Z\ + 
+ (12a^6 + 8a6^ - 6ab) Zf ^ + 4a6^ Z\ + 406^ Z^^\ + 
+ [60^62^2 + 6a V]Z2-2, 



(150) 



where the Z™ at the right-hand side of Eq. (150) should be evaluated at 



p'e^'''. 



7.1.2 Computation for Z\ 

We have 



z\{a^b(^e'>') = Y^ K\^:z:::{p'^'') 



(151) 
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where we have K^!^ Z™ ^ only if n' and m' have the same parity and 
|r?7,'| < n' < 3 — |r7i' — 1|. This leaves us with the cases 

m' = —1, n = 1 ; m = 0, n' = 0,2 ; m = 1, n = 1,3 ; m = 2 , n = 2 . 



Thus we compute subsequently 



(152) 



m' = —1. a. K^[ = T^l — T33 by Eq. (26) with T^l given by Eq. (27) 
upper case, as 



^ilP =l,q = 2;p' = l,q' = 0] = ?,a% P','^'\l - 2A^) P','^'\2B' - 1) , 



2ap(2,1), 



,2n p(2,1)/or2 



and T33 = since n — n" = < 4 = |m — m"\, see Eq. (27) 



(153) 



m' = 0. a. Kl^ = T^^ - T^^ with T^^ and T.^^ given by Eq. (27), upper 



CcLSG, cLS 



(1,0), 



i2n d(1>0)/od2 



r3^o> = l,g = 2; p' = 0,g' = 0] = aPr''^^(l - 2A') Pr\2B' - 1) , 



2 d(1,2). 



l2\ d(1.2)/od2 



(154) 

Tl!\p = 1, g = 2 ; p' = 1, g' = 1] = 3a6^ Po^^''^(l " 2A^) Po^'''^(25' - 1) • 

(155) 

m' = 0. b. ^3^^ = T^^ - T^^ with T^^ given in Eq. (154) and T^^ = since 



n — n" = —I < I = \m — m"\, see Eq. (27). 



m' = 1. a. K^l = T3Y — T33 with T^^ and Tgg^ given by Eq. (27), both upper 
and middle case, as 



T^llp = l,q = 2;p' = 0,q' = l] = bPr'{l - 2A')Pr'{2B' - 1) , 



(0,1) 



\2^ p(0,l)/r)R2 



(156) 



T^^lp = l,q = 2;p' = l,q' = 2] = b^ P^'il - 2A') Pr'{2B' - 1) . 



,3 p(0,3). 



i2n p(0,3)/r)R2 



(157) 



m 



' = 1- b. Kll = T3V - T3i5^ with T3I3I given by Eq. (156) and T^^ = 0. 



m' = 2. a. i^gl = T^2 ~ '^il with T^2 given by Eq. (27), middle case, as 
T^^[p = l,q = 2;p' = 0,q' = 2] = 3ah' Pi,'''\l-2A') P^'^'\2B'-1) , (158) 



and Til = 0. 
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We calculate the right-hand sides of Eqs. (152)-(157) as 



This gives 





= 3a% , 






(159) 


rplQ 

-'so 


= a{3a^ + Qb^ - 2) , 


rplO 


= Safe^ , 


(160) 


rpll 


= b{Qa^ + 3b^ - 2) , 


rpll 


= h\ 


(161) 


rpl2 
-'32 


= 3ab^ . 






(162) 


Klf 


' = 3a% , 










= 3a^ + 3ab'^ - 2a , 


1^10 . 


= Sab"^ , 


(163) 


-f^ai 


= 6a% + 2b^ - 26 , 




= &^ 


(164) 


1^12 

-'^32 


= 3a&2 . 









Hence 



Zl{a + bp' e'"^') = 3a% Zf ^ + [(Sa^ + 3ab^ - 2a) Z° + 3ab^ Z°] + 



+ [{6a% + 2b^ - 2b) Z\ + \? Zl] + 3ab'^ Z 



2 ry2 

2 ' 



(165) 



where the Z™ at the right-hand side of Eq. (166) should be evaluated at 



p'e^'''. 



It is obvious that for the expansion of a general ^^(a + b p' e**' ) one can 
construct a concise and efficient computer code on basis of Theorem 2.1, 
taking advantage of the various shortcuts and reuse of intermediate results 
such as those encountered in passing in the above two examples. 

7.2 Example for Section [3] 

We compute, using Theorem 3.1 and the notes thereafter, the Zernike 
expansion of Z^ **corr Zq. It is easy to show by elementary means, Z^ **corr Zq 
being the area of the common part of two disks of radius 1 whose centers are 
at a distance p apart, that 



\,Zq * *corr 



Z°)(p) = 2 arccos(|p)-^pJl-(ip) 



< p < 2 . (166) 
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From Eq. (62) we have 



oo 
(Z° **corr^o)(p)=27r / 



Jfju) Jo{pu) 



u 



du . 



(167) 



u 

This integral can be found in [27], 2.12.42, item 31 on p. 232, and this would 
yield Eq. (166) when the parentheses would have been placed correctly in 
this reference (a cross-check with [37], 2.12.42, item 15 on p. 230, that arises 
when the integral on the right-hand side of Eq. (167) is differentiated with 
respect to p, shows inconsistency of [37j in this matter). 

We have m" = m — m' = and n" is even at the right-hand side of 
Eq. (63), and this yields 

(^O°**con-^0°)(P) = ^ E (^"+l)rZ"^°"(ip), 0<p<2 



fl68) 



n" even,>0 



where 



oo 

rX" = 8(-l)-^""|^iM^iH^n"+i(2n)^ = 

n 



= 2(-l 



)i""[QSo+^ + 2QS^^ + Q^;'+^] 



(169) 



with 






n" = 0,2,... , (170) 



n" = 2,4,... , (171) 



rrl„/M1^2 32^i(n"-2)' n -^,<0,... , {III} 



^" ((in")! 



while (5o2 ^^d Q\2i Q22 vanish. Thus this yields the Zernike ^-expansion of 
the function 

J2 Cl„ Zl„{T = \p) , 0<r< 1.(173) 

fivpn.>n 



2(arccos t — r a/1 — t' 



n" even,>0 

The C°„ are given in integral form as 



1 

Cl„ = 4(n" + 1) /" (arccos t-t VI-t^) RI„{t) t dr ; (174) 
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the evaluation of the integrals in Eq. (175) becomes cumbersome, already 



for low values of n" = 0,2,.... We compute from Eqs. (169)-(171) and 



r^O _ TT / // I 1 N -nooo 
^n" — aV'' ^ ^) ^ OOn" 



TT 






8 ' ^ 32' 



(175) 



Unfortunately, there does not seem to exist a closed formula for the values of 
-Pfc ' (0) ^s required in Eqs. (169)-(171), except for the case 7 = 5 = 0, see 



Table 22.4.1 on p. 777. Furthermore, the C°„ decay only slowly because 
of non-smooth behaviour of (Zq * *corr ^o)(^)/^) around {y^jj) = (0,0) and, 
to a lesser extent, around z/^ + /x^ = 1. 



7.3 Examples for Section [4] 



We have computed, using Theorem 4.1, the Fourier coefficients a^ in the 
cosine representation 



Ln/2J 

R^{cosx) = y^ cin-2j cos{n — 2j) X 

3=0 



(176) 



for various cases of integer n, m with n, m non-negative and n — m even and 
non-negative. The results are collected in Table I. 





(a) 



(b) 



Figure 6: Plot of (a) RIq{p),() < p <l, and (b) R%{cosx), < a; < 7r/2. The 
sampling rate used to display (a) is not high enough to adequately represent 
the last peak but one just before p = 1. 



TABLE I. 
R^ in polynomial and cosine representation with p = cos x and Ck = cos kx 
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i?4 
Rl 

Ra 



Ry. l = co 

R9, : 2p2-l = c2 

i?o : 6p4-6p2 + l = lco + |c4 

i?o : 20p6 - 30p4 + 12p2 _ 1 = | C2 + | Cg 

i?o : 70p8 - 140p6 + 90p4 - 20p2 + 1 = Jl Co + ^ C4 + I C8 

i^i : p = Ci 

i?i : 3p3-2p=ici + |c3 

i?i : 10p5-12p3 + 3p=|ci + |c3 + |c5 

R'r : 35p7-60p5 + 30p3-4p=J|ci + |c3 + ^C5 + ic7 

-^2 • P = 2 '^0 "*" 2 '^2 

4p4 _ 3p2 = 1 C2 + I C4 

15p6 - 20p^ + 6p2 = ^ Co + ^ C2 + ^ C4 + II C6 

56pS - 105p^ + 60p^ - 10p2 = ^ C2 + j^ C4 + ^ C6 + ^ C8 

Rl : 5p5-4p3 = |ci + ^C3 + ^C5 

i?3 : 21p7-30p5 + 10p3 = lfci + ^C3 + ic5 + ic7 

-^4 • P = g Co + 2 ^2 + 8 C4 

i?| : 6p6-5p4 = ^C2 + ^C4+^C6 

i?^ : 28pS - 42p6 + 15p^ = ^ cq + 3^ C2 + | C4 + ^ Ce + ^ cg 

Rl : P' = ici + ^C3 + ^C5 

Rl : 7p^-6p5 = Aci + ||c3 + ic5 + ^C7 

-^6 • P ~ Te '^0 + 32 C2 + Jq C4 + 32 ^6 

_Rg : 8p — 7p = 32 '^2 + Yg C4 + 32 ce + xg cg 

-^7 • P = 64 '^1 + 64 '^3 + 64 '^5 + 64 '^7 

-^8 • P = 128 '^0 + 16 '^2 + 32 C4 + 16 C6 + Y28 ^8 

In Fig. 6 we have displayed i?gQ(p), < p < 1, and -RgQ(cosa;), < x < 7r/2 
to illustrate the point that the variation of the radial polynomial is more or 
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less spread out uniformly over the x-interval. 



A Convergence of the integral in Eq. (7) 

We shall show in this appendix that for non-negative integers n and m 
with n — m even and non-negative, the integral 

V 

lim / Jn+i{u) Jmipu) du (177) 



11— >oo 





converges to (—1)^^ Rn{p) for < p < 1 and to for p > 1, and that it 
does so boundedly in p > and uniformly in p outside {1 — e,l + e) for any 
£ > 0. We have from il4|, 9.2.1 on p. 364 

/ 2 1 1 

Jj.[u) = \ — cos(m- -A;7r- -tt) + 0(m~^/^) , m -^ oo , (178) 

V VTM 2 4 

and Jkiu) is smooth and bounded on m > 0. Therefore, to show bounded 



and uniform convergence of the integral in Eq. (177) on the appropriate sets 



of p, it is sufficient to establish this for the integral 

2 r 1 , 1, X 1 X / 1 1 X , / X 

lim — / — cos(m — -(n + 1) TT — -tt) cos(pM — -rriTT — -tt) ctM . (179) 
v^co 71 J u 2 4 2 4 







Once this has been established, the issue of to what the integral in Eq. ( 177) 
converges is settled by the remark that Z,^ and 27rz" exp{irrnp) J„+i(27rr)/27rr, 
see Eq. (Is]) , are 2D Fourier pairs so that Fourier inversion of the latter func- 
tion yields the former in L^(M)-sense while the former function is smooth 
outside the set u'^ + fj? = 1. 

Using elementary trigonometric identities, we have 

2 ? 1 / 1/ .N In/ 1 1 N J 

— / — COS n [n + ljiT TTjcosipu rmr ttjcLu = 

TT J u ^ 2^ ^ A ^ ^^ 2 4^ 

1 

■1)P ?sin(l-p)M, (-1)'' } cos(l + p)u, 

du / du = 



JC ,1 U TC J u 

1 



{l^p)v {l+p)v 

'')' f ^^"^dx-tll! f 'J^dx, (180) 



TT J X IT J X 

1-p 1+p 
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where we have set p = ^{n — m), q = ^ {n + m). Since both functions 

y y 
sin a; , f cosx , 
dx, y>0; / dx , y>l, 181 

X J X 

1 

are bounded and have a finite hmit as y — )■ oo, the convergence of the integral 



in Eq. (179) is bounded in p > and uniform in any closed set of p's not 
containing 1. The assumption that n and m have same parity is essential: 

oo 

f - cos^ u du = OO. 

J u 

1 
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